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Nuclei in dense matter are influenced by the medium. In the cluster mean field approximation, 
an effective Schrodinger equation for the ^-particle cluster is obtained accounting for the effects of 
the surrounding medium, such as self-energy and Pauli blocking. Similar to the single-baryon states 
(free neutrons and protons), the light elements (2 < A < 4, internal quantum state v) are treated 
as quasiparticles with energies EA,v(P',T,n n ,n p ) that depend on the center of mass momentum 
P, the temperature T, and the total densities n„,n p of neutrons and protons, respectively. We 
consider the composition and thermodynamic properties of nuclear matter at low densities. At low 
temperatures, quartetting is expected to occur. Consequences for different physical properties of 
nuclear matter and finite nuclei are discussed. 

PACS numbers: 21.65.Ef, 05.70.Ce, 25.70.-q, 26.60.Kp, 26.50.+X 

INTRODUCTION 

Quantum condensates are one of the amazing phenomena in many-particle physics. Isospin-triplet (proton-proton 
or neutron-neutron) pairing is well established [l| not only in nuclear matter at low temperatures, but also in finite 
nuclei. Isospin-singlet (neutron-proton) pairing may become of relevance in symmetric and asymmetric nuclear matter 
0. On the one hand, the interaction is stronger, even a bound state, the deuteron, can be formed. This would lead 
to higher transition temperatures and the transition from Cooper pairing to Bose-Einstein condensation. On the 
other hand, with increasing difference of the chemical potentials and increasing Coulomb effects, the formation of a 
condensate becomes more difficult, in particular in asymmetric matter. The strong interaction in the isospin-singlet 
channel leads to a precursor for the Bose condensate, the occurrence of a pseudogap [3j above the critical temperature. 
Furthermore, pairing competes with the formation of an a-particle condensate (quartetting) in particular at low 
densities. Therefore, the role of correlations, in particular the formation of clusters and of a quantum condensate, is 
an interesting item in strongly interacting quantum liquids such as nuclear matter. 

Systematic quantum statistical approaches, see have to be used for a treatment using appropriate concepts 
such as thermodynamic Green functions, spectral function, and self-energy. Within a chemical picture, a cluster 
decomposition of the self-energy can be performed that considers the formation of clusters in a dense medium, 
taking into account symmetrization (Pauli blocking), screening of the interaction, dynamical self-energy, and other 
many-particle effects. An interesting approximation is the cluster mean- field approximation 0-0] that considers the 
few-particle Schrodinger equation in a correlated medium, where the correlations in the medium have to be determined 
in a self-consistent way. This would allow to describe nuclear matter in a wide range of densities, from saturation 
density where a mean-field approach on the single nucleon level is possible, to the low-density region where the nuclear 
statistical equilibrium (NSE) applies and a-matter can occur. We discuss the formation of clusters in nuclear matter 
in Sec. 2 and give some applications to finite nuclei in Sec. 3. 

NUCLEI IN MATTER 

The few-body problem describing A < 4 nucleons in hot and dense matter can be related to an in-medium wave 
equation (Bethe-Salpeter equation) that is derived from many-particle approaches, see Ref. We consider only 
bound states v, dropping the spin quantum number. For the light elements, v — d,t,h,a denotes the deuteron ( 2 H), 
the triton ( 3 H), the helion ( 3 He), and the a particle ( 4 He). Considering uncorrelated nucleons in the medium, the 
few-nucleon wave function and the corresponding eigenvalues follow from solving the in-medium Schrodinger equation 

[Ef(l) + ... + Ef(A)]^ P {l...A) 
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For brevity, the single- nucleon quasiparticle energy E^(pi) is denoted as E^ p (l). The nucleon- nucleon interaction 
V{ij,i'j') becomes medium dependent due to the Pauli blocking prefactor [1 — fi(i) ~ /i(i)]- The phase space 
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TABLE I: Parameter values for the Pauli blocking shift A£^ auli (P; T, n, Y p ), Eq. ©, in the low-density limit 
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occupation is described by a Fermi distribution function normalized to the total density of nucleons, 
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in the low-density, non-degenerate limit (p T < 0). The chemical potential /i r is determined by the normalization 
condition 2f2 _1 fi(p) = n T , where r denotes isospin (neutron or proton), and has to be expressed in terms of these 
densities n n , n p or the baryon density n = n n + n p and proton fraction Y p = n p /n, and the temperature T. 

The in-medium Schrodinger equation ([1]) contains the effects of the medium in the single nucleon quasiparticle shift 
as well as in the Pauli blocking terms. Obviously, the bound state wave functions and energy eigenvalues as well as 
the scattering phase shifts depend on temperature and density. In particular, we obtain the cluster quasiparticle shifts 



E?(P; T, n, Y p ) - E U {P) = AE^(P; T, n, Y p ) + AE^ ull (P; T, n, Y p ). 



(3) 



The contribution of the single nucleon energy shift to the cluster self-energy shift A_E^ E is easily calculated in the 
effective mass approximation, where the single-nucleon quasiparticle energy shift A£f E (l) can be represented by the 
energy shift AE^ E (T, n, Y p ) and the effective mass m*(T,n,Y p ), see 0. Different expressions for the single-nucleon 
quasiparticle shifts are available such as Skyrme, RMF, or DBHF, see Q. Since the single-nucleon quasiparticle 
energy shift arises in the continuum of scattering states as well, the influence on the binding energies is small. 

We consider here the Pauli blocking shift of the binding energies 



AE^ U[1 (P; T, n, Y p ) = E™(P) - E V {P) -{A- Z)AE™ - ZAE 



(4) 



Different fit formula have been given [H, 0, H3] to parametrize the solution of the in-medium Schrodinger equation (JTJ . 
We follow the recent work [8|. 

In the low-density limit, a linear dependence of the energy shifts on the nucleon density n follows from perturbation 
theory. For a more general dependence on the total nucleon density we consider the expression 



AE^\P- n, T, Y p ) = c v {P; T) < 1 - exp 



Cv{P;T) 



y v (Y p )n - d v (P;T, n)n' 
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where the dependence on the asymmetry is given by yd{Y p ) — y a (Y p ) = 1, for triton yt(Y P ) = (§ — f^p)> an< 4 f° r 
helion yh{Xp) = (f + f^p) * na * renec ts the different proton and neutron content.. The functions /„(P;T,0) can be 
calculated in first order perturbation theory using the unperturbed wave functions of the free nuclei. We performed 
model calculations with a separable potential that reproduces the empirical binding energies and rms radii. Motivated 
by the exact solution for A = 2, we use the following fit for arbitrary v: 



f„(P;T,n) = /„,iexp 
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For zero momenta, P = 0, the temperature dependence of c„(0; T) and d v (0;T,n) is expressed as 



c u (0;T) = c Ul a + 



(T - c„ :2 ) 2 + c„3 ' 



d„(0; T, n) 



(7) 
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TABLE II: Parameter values for the Pauli blocking shift A£^ auli (0; T, n, Y p ), Eq. ©, in units of MeV and fm (c„, , c„, 2 , d„, 2 - 
[MeV]; c„,i - [MeV 3 ]; c„, 3 ,d„,3 - [MeV 2 ]; d vA - [MeV 2 fm 6 ]; u v - fm; v u - [MeV" 1 fm]) 
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The Pauli shift at finite momenta is fitted with c !/ (P;T) = c^(0;T) not depending on P, but 

d v (P;T,n) = d v {0;T,n)e-y& . (8) 

An additional dependence on n, T is considered at finite values of P. Another additional dependence on n, T for finite 
momenta is introduced in /„(P;T, ri) where the dispersion relation becomes density dependent due to the parameter 
tt„. Parameter values are given in Tab. fl] and Tab. [TT1 

The in-medium Schrodinger equations ([1]) describe in the case A = 2 (d) and A = 4 (a) also the onset of a quantum 
condensate if according to the Thouless condition the in-medium binding energy coincides with Ap. 

With these parameter values, it is possible to evaluate the composition of warm nuclear matter at arbitrary baryon 
density n below the saturation density, temperature T, and proton fraction Y p . Equations of state can be obtained 
for different thermodynamic properties, see also Q. A smooth transition from the description of nuclear matter near 
saturation density to the low density region where the NSE can be applied has been found. A systematic inclusion of 
mean-field effects is possible. 

An unsolved problem is the behavior of the internal energy at low temperatures. The internal energy for symmetric 
matter is shown in Fig. 1. At low densities and temperatures it approaches the binding energy of the a particle. 
With increasing density, the self-energy shift leads to the behavior dU/dn < at T — that means thermodynamical 
instability. With the disappearance of the bound states, the internal energy approaches the values for the quasiparticle 
single-nucleon approximation. This occurs rather abruptly near n = 0.02 fm~ 3 so that an intermediate region 
of metastability appears. However, this may be changed if better approximations for the internal energy in the 
intermediate region are available. Alternative approaches to the internal energy of symmetric matter at T = have 



been considered in connection with the formation of a quantum condensate [111 . [1 



a CLUSTERING IN NUCLEI 



Nuclear matter at low densities and temperatures is thermodynamically instable. This region of phase instability 
can only be reached in nonequilibrium or in inhomogeneous systems. An interesting application are nuclei, where 
nuclear systems at low densities can occur. In particular, excited states of symmetric 4n nuclei near the na threshold 
may show an a cluster state such as in 8 Be, 12 C, possibly 16 0, and others. As a simple ansatz to describe a clusters 
with the same c.o.m. orbits, like the pairing state, the THSR wave function was introduced This wave function 
is a product ansatz for the c.o.m. motion of antisymmetrized a clusters, similar to the BCS wave function. This 
approach has been proven to describe different properties of such low-density states successfully. 

In particular, a relation between the density and the condensate fraction has been discussed [14| that indicates the 
disappearence of the condensate around n sa t U ration/3. Whereas a clusters are well developed in the Hoyle state, in 
the ground state of 12 C no a clustering is present and a quasiparticle single-nucleon shell model is applicable. The 
reason is the antisymmetrization of the nucleonic wave functions that results in the Pauli blocking. In contrast to the 
calculations shown in Fig. 1, where the Pauli blocking is caused by uncorrelated nucleons that form a Fermi sphere, 
the Pauli blocking by a clustered medium is less efficient because the nucleonic wave function is more distributed in 
momentum space. Therefore, we expect that the crossover from a matter behavior, given by a cluster- virial expansion, 
to the quasiparticle nucleonic liquids happens not near p — 0.02 fm -3 , but near p = 0.05 fm -3 so that a continuous 
decrease is expected. Then, the intermediate region of metastability where dp,/dn > would disappear. 
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FIG. 1: Internal energy U for symmetric nuclear matter at different temperatures. 



The behavior of the quasiparticle shift E% p (P;T,n,Y p ) — E a (P), Eq. is sensitive to the values of self-energy 
shift and Pauli blocking that partially compensate. At zero temperature, the self-energy shift amounts about -4000 n 
MeV fm 3 , whereas the Pauli blocking term is larger than 4434 n MeV fm 3 as given for T = 1 MeV in [5j. If we estimate 
the kinetic energy of the nucleons bound in the a particle as E^ m = 51. 6 MeV [8| and simulate the corresponding 
occupation in momentum space by an effective temperature E^ ln /6, the Pauli blocking shift is reduced by a factor 2 
so that the total quasiparticle energy shift becomes negative, around -2000 n MeV fm 3 at T = 0. A decrease of the 
internal energy follows also from the cluster virial expansion that gives about -1500 n MeV fm 3 at T = 1 MeV. 
This value is only based on empirical a — a scattering phase shits. The Ali-Bodmer interaction, which is also based 
on these phase shifts, has been considered in [12J to calculate the a condensate fraction at T = 0. Also a decrease 
of internal energy was found. Exploratory lattice calculations for a T = a condensate state [ll| gave no definite 
answer. Whereas for 4 He an increase of energy with increasing density was obtained, the opposite occurs for 16 0. 

Another interesting region where a clustering is expected to occur is the surface of nuclei. Like pairing that is 
increasing near the surface [l[ where the density drops down from the saturation value, also a clustering and possibly 
quartetting can occur. A simple way to incorporate this is a local density approach. Estimations of the role of isospin- 
singlet pairing and quartetting for nuclei near the N = Z line have been given in [l6| analysing binding energies of 
the ground state and of the excited state. The Wigner energy can be explained this way. This effect is seen in the 
light nuclei where the shift of the proton chemical potential due to the Coulomb energy is not strong compared to the 
condensation energy. 

Cluster formation at the surface of nuclei is of interest also in the region of superheavy elements. As well known, 
the most dominant decay channel is a decay, in contrast to the emission of protons or neutrons. Therefore one expects 
that a particles are preformed in the skin of the superheavies where the density drops down. A corresponding picture 
is a stable core such as the double magic lead or uranium nucleus, surrounded by a skin of nucleons at reduced 
density where cluster formation is possible. Of course, the Pauli blocking mechanism has to be considered beyond the 
local density approximation, considering the density matrix due to the core orbitals. We discuss two consequences 
of clustering in the skin of superheavies, the binding energy and the rms radii. Further signatures can be seen in 
reaction processes. 

The staggering of the binding energies of superheavies is well known from experimental data. It would be of 
interest to what extent this staggering can be explained within a quasiparticle single-nucleon approach such as shell 
calculations, or using a THSR-like wave function that describes the c.o.m. motion of antisymmetrized clusters. 

Another effect are the rms radii of superheavy elements. Experimentally, the charge radii (r 2 ) c harge are measured. 
The calculated values are the point radii (r 2 ) po int that consider point-like nucleons. The relation between both is 
given by the relation 17 1 



(r ) P oint + (r ) proton + (N/Z)(r 2 ) neutron — ) charge 



(9) 
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that takes into account the contribution of the proton form factor due to its spatial extension, (r 2 ) pr oton = 0.743 fm 2 , 
and the neutron contribution (r 2 ) nout ron = —0.116 fm 2 . If we assume that the surface is formed by a-like clusters, we 
expect another relation, 

(?" 2 }point + (r 2 )a = (r 2 ) ch argc (10) 

with the charge rms radius for the a particle, {r 2 ) a = 2.822 fm 2 . Assuming that the point rms radii are nearly the 
same, the effect of a-clustering would increase the rms radius according (r 2 ) a — (r 2 ) plot0 n — {N/Z)(r 2 ) ncutr0 n = 2.195 
fm 2 . For instance, considering a heavy nucleus with [(r 2 ) point] 1 / 2 = 6 fm, the charge rms radius would be larger 
by 2.98 % if the a cluster determine the surface, compared to the usual nucleonic surface. The empirical values of 
charge rms radii of nuclei [13] have been compared with Hartree-Fock Bogoliubov shell model calculations [HJ], and 
excellent agreement has been found with exception of the superheavies like Cm where the measured values of the 
charge rms radii are about 2 % higher than the calculated ones. An interesting point would be whether the difference 
can be explained by clustering. Assuming this, a emitters should have an expanded charge rms radius, compared to 
quasiparticle single-nucleon calculations. 
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